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Abstract 

Nonlinear waves in a liquid containing gas bubbles are considered 
in the three-dimensional case. Nonlinear evolution equation is given 
for description of long nonlinear pressure waves. It is shown that in 
the general case the equation is not integrable. Some exact solutions 
for the nonlinear evolution equation are presented. Application of the 
Hirota method is illustrated for finding multi-soliton solutions for the 
nonintegrable evolution equation in the three-dimensional case. The 
stability of the one-dimensional solitary waves is investigated. It is 
shown that the one-dimensional solitary waves are stable to transverse 
perturbations. 



1 Introduction. 

In the nature, technology, and medicine we often can observe the medium 
consisting of a liquid with gas bubbles [l}j3]. One of the important problem 
for investigation in bubbly liquids is propagation of waves. Wave processes in 
bubble-liquid mixtures were considered in many papers. Taking into account 
the linear approach we have the well known problem that does not differ in 
essence from many other problems of the linear wave mechanics. 

However a lot of new possibilities appear when we consider the nonlinear 
processes. The influence of nonlinearity on waves propagation in a liquid 
with gas bubbles was first studied in |4|j6]. It was shown that long weakly 
nonlinear waves in a bubbly liquid can be described by the Korteweg-de 
Vries equation. Authors also took into consideration the viscosity of liquid 
and obtained the Korteweg-de Vries-Burgers equation for the description of 
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nonlinear waves. However in [4j-[6] only one-dimensional weakly nonlinear 
waves were considered. 

Recently [7J|8) we have studied the influence of heat-transfer between a 
liquid and a gas in bubbles on propagation of weakly nonlinear waves in 
bubble-liquid mixtures. As a result we have obtained nonlinear evolution 
equations of the second, third and fourth order for the description of wave 
processes in a one-dimensional bubble-liquid mixture. 

The real nonlinear waves in a liquid with gas bubbles are two- and three- 
dimensional. So we need to use multidimensional mathematical models for 
studying nonlinear waves in a bubbly liquid. Some problems of pressure 
wave distributions in a bubble-liquid mixture in the two-dimensional case 
were considered numerically in (9 11 



However the calculations of wave processes presented in these works are 
based on cumbersome mathematical models and do not allow us to reveal 
basic parameters of nonlinear wave processes in a bubbly liquid. Besides we 
cannot find any exact solutions for these models and consequently cannot 
effectively analyze correctness of numerical modeling. 

Thus it is important to consider long weakly nonlinear waves in a liquid 
with gas bubbles in the three-dimensional case. In this case we can obtain 
evolution equation for the nonlinear waves in a bubbly liquid. The evolution 
equation can be studied both analytically and numerically. Also analyzing 
the evolution equation we can determine nature of physical characteristics 
influence on nonlinear waves. 

The aim of this work is to study long weakly nonlinear waves in a liquid 
containing gas bubbles in the three-dimensional case. In section 2 we present 
governed system of equations for nonlinear waves in a bubbly liquid. In 
section 3 with the help of the reductive perturbation method we obtain the 
evolution equation for long weakly nonlinear weakly three-dimensional waves 
in a liquid with gas bubbles. In section 4 we show that this equation is not 
integrable in the general case. In section 5 we obtain some exact solutions 
for the evolution equation. We present stability analysis of one-dimensional 
solitary waves to transverse perturbations in section 5 as well. 



2 System of equations for the description of 
wave processes in liquid with gas bubbles. 

Let us suppose that a liquid containing gas bubbles is a homogeneous medium, 
which has an average pressure and a temperature [TJ[2] . We assume that the 
liquid is incompressible and bubbles are spherical. We do not take into con- 
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sideration formation, destruction, interaction and coalescence of bubbles. Let 
us assume that all gas bubbles have the same size and the amount of bub- 
bles in the mass unit is constant N. We suppose that the gas in bubbles is 
ideal and pressure of gas in bubbles is governed by polytrophic law. We shall 
consider the viscosity of the liquid only in the interphase boundary. Under 
these assumptions the following closed system of equations may be used for 
the description of waves processes in the liquid containing gas bubbles [TJ[2] 

^ + divpu = 0, (1) 



p ( ^ + uVu 1 + VP = 0, (2) 



Pl [ RR TT + |i£ + ^R T \ = Pg -P- 2 ^ : (3) 



p \ 3n 

H 



P 9 = P 9 ,o[^) > (4) 



P= i 1 ~ <t>) Pl + <t> Pg, (5) 



<P = Vp, V=-ttR 3 N (6) 
3 

where £, r], ( are cartesian coordinates, rj, (, r) is the density of the bubble- 
liquid mixture, P(£, 7/,£,t) is the pressure of the mixture, u(£, 77, £, r) = 
i], C, t), u^(^, rj, (, r), u^ 3 \^, 77, (, r)) T is the velocity of the mixture, 
R = P(£, 77, C> t) is bubbles radius, pi, p g (£, rj, (, r) are densities of the liquid 
and the gas respectively, P s (£, 77, £,r) is the pressure of the gas, P 9i0 and R 
are the pressure of the gas and the radius of bubbles in the unperturbed 
state, a is the surface tension, v is the kinematic viscosity of the liquid, 
is the volume gas content, V is the gas volume in the unit mass of the 
mixture, n is the polytrophic exponent and the nabla operator is given by 

v = (d/d^d/d v ,d/dC). 

We shall assume that in the unperturbed state the pressure and the den- 
sity of the bubble-liquid mixture are constants, bubbles have all the same 
radius and uniformly distributed in the liquid. 
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Assuming that the mean inter-bubble distance d is much greater than the 
bubble radius, we do not take into account bubbles interaction. We suppose 
that the characteristic wavelength L is greater than linear scale of a region 
containing many bubbles [IJ[2j[4j[6]. Thus we have the following relation 
between characteristic length scales in the bubbly liquid 

Ro < d < L (7) 

Assuming that the volume gas content is small C 1 we have from ([5j 
and @ 

Pi fo \ 



' 1 + PiV 
Using (|8]) system ([T])-(|6]) can be written as 

^ + divpu = 0, (9) 



J^ + uVu) +VP = 0, (10) 



P= T - A 17 , V=Kr 3 N. (12) 
1 + Pz V 3 

Let us suppose that deviation of the density from the unperturbed density 
is small in comparison with the last 

p(£,ViC,t) = po + p(£,ri,C,T), p = const, \\p\\ « p , (13) 

where po is the density of the bubble-liquid mixture in the unperturbed state. 

Taking into account (13) from (12) with the accuracy up to terms of order 
p 2 we have 

R = R - pp + pip 2 , 

p3 _ 3 (I 1_\ Rp R (3poV -l) 

°"4vrivU pj' ^'SplVo' ^~ 9p^ 2 ' (14) 

Vo = UNI& 



4 



Substituting (14) into (J9j) ( 12 ) and taking into account (13) we obtain 
dp 



dr 

(Po + p) 



+ podivu + divpu = 0, 

uVu] + VP = 0, 



Pn 



du 

Ih 

2a ( 3n \l P 2a p 

P V R 



F? 2 



P+ 



+ 



3n[(3n + 1) p 2 - 2 ^ P ]^o c^(p 2 - Pi^o] 



(15) 



+Pl p PoPrr - Pi (2pi P + P )P Prr - Pi ( 2p X P, 

4z/pp, _ (4z/p 2 - 8up 1 R )pi _ „ 



p 2 + 



3P 2 
2 



^2 



3Po 



Pr + 



3-Ro 



PPt 



where we denote by Po the pressure of gas in bubbles in the unperturbed 
state. 

Linearizing system (15) and assuming that P is proportional to p we 
obtain the following equation 



c Ap, c 



2 3npP 2aji 



Pn 



d2 ' 

K 



(16) 



where A is the Laplace operator. From ( 16 ) we can see that speed of waves 
in the linear case is cq. Further we shall use Co for introducing dimensionless 
variables. 

Let us introduce the dimensionless variables 



£ = L£, V = L V , ( = L( r 



L 

Co 



T 



U = Cq u 



P = Pop , 



P = PnP / + Pn- 



2a 
Po 



(17) 



where Po— 2a /Rq is the pressure of the mixture in the unperturbed state. The 
quantities L and r* = L/cq are characteristic length scale and characteristic 
time on which parameters of the bubbly liquid change 12]. 

Using the dimensionless variables we can reduce (15) to the following 
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system of equations (the primes are omitted) 
dp 



dr 



+ divu + divp u = 0, 



+ ( <^ + uVu 



1 



VP = 0, 



a 



;is) 



ap + ocxp + (3p TT - (pi + /3 2 )pp T T - (Pi + -/3 2 )/V+ 



+xp r + xipp r 



a 



3npp 2<rppo 



Rq 



RqPq 



ft 



Pi pRpclpo 



a i 



3n[(3n + 1) p 2 - 2 p x R ]p 2 a(p 2 - p 1 R )p 2 



2R 2 



PoR 3 



Pi 



2pi piR c 2 Q p 2 



P L 2 



Pi 



PnL 2 ' 



(19) 



Aup pipoco 



Xi 



(4z//i 2 - 8vp 1 R )p l p 2 l c 



3R 2 P n L 



3RqPqL ojxqI'o^ 
Further we shall analyze this system of equations with the help of the 



reductive perturbation method. From (18) we shall obtain the nonlinear 



evolution equation for long weakly nonlinear weakly three-dimensional waves 
in the liquid with gas bubbles. 



Nonlinear Evolution Equation for long waves 
in the liquid with gas bubbles in the three- 
dimensional case. 



Let us consider transverse variations of the long waves governed by (18). We 



suppose that waves propagate along the with allowance to transverse 

evolution in y and z directions. To obtain the nonlinear evolution equation 
for the long waves we shall use the reductive perturbation method, which was 
introduced in 12,13 . The application of the reductive perturbation method 
to multidimensional problems, for example, were considered in 14 17 . For 
the modern review on the reductive perturbation method see 18 



In accordance with reductive perturbation method we need to introduce 
scale transformations of the independent variables. To this end we use the 
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dispersion relation for the linearized system (18), which in the case of x = 
is the following 



kg ~\~ k^ -\- k^ 
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(20) 



where kg, k v and kg are wave numbers and a; is a frequency. 

As we consider weakly three-dimensional waves we neglect wave numbers 
k„ and kg in (20). Assuming that w 2 < 1 form (20) we obtain 



P 3 

k e ~ u + — oj 
2a 



(21) 



Assuming that u 2 = ep 2 we see that phase of plane wave in the x-direction 
has the form 



kgi-ur=~^e^{i-r) + ^e^ 



Using (|22j) we can introduce 'slow' variables in (|18|) as follows 

x = e"\i - r) 



t 



,3/2 



r, y = eSr], z = e^C > £< 1. 



(22) 



(23) 



Then we have the following relations 







d£ dx dr dt dx dr) dy d( dz 



In (23) we suppose that transverse variables y and z are 'slower' than x 



as far as we consider transverse perturbation of the wave propagating in the 
x direction. 

We consider the long waves in the liquid with gas bubbles. Thus the char- 
acteristic wave length is much greater than the mean inter-bubble distance 
and the radius of bubbles in the unperturbed state (see (|7|). In this case we 
can interpret small parameter e as e — <C 1. Also we assume that 5 and 
7 are parameters which characterize perturbations in y and z directions. 

Substituting (23) in (18) and dividing on e 1//2 in first three equations we 



have the system of equations in the form 

ept - p x + u ( x 
+e l 'H~ P u^) z = 



e~ Pt - ~ Px + + e^buf + e^Tuf ) + + e 1 / 2 5(pu^) y + 



(25) 



1 + p){eu? - v$ + vg> + e^Su® + e x '*>yvP> 




a 



(26) 



7 



(1 + p)(eu {2) - uf + u« uf + e l ' 2 5u^ wf + e^ya® u®)+ 

i/ 2 5 ' ( 27 ) 

+£ 1 / 2 -P y = 

a 



(1 + p)( £ uf } - uf + u« uf + e l ' 2 5u^ + eVV 35 uf)+ 

i/ 2 7 ( 28 ) 



P = ap + a x p +e A Pptt + E (3 Ptx + £ (3 Pxx ~ £ (Pi + ^)pptt+ 



2(/3i + /3 2 )pp te + #0/5/5** - £ J (/5i + -/3 2 )p7+ 

3 3 
+e 2 Wi + ^ 2 )pxPt - e (ft + -/3 2 )p 2 + £ 3//2 xpi - £ 1/2 xp,+ 

+e 3/2 X!ppi -e 1/2 xpp, 



(29) 



Let us look for the solution of (|25|)-(|29|) in the form 



u x ' = eu 

(3) 



u (2) = £ 3/2 M (2) + £ 5/2 M (2) + 



£ 3/2 M (3) + £ 5/2 u (3) + 



p = epi + e 2 p 2 



(30) 



P = eP 1 + e 2 P 2 + 



Substituting ( |30| ) into (J25J)— (|29j) and equating expressions at e we have 

(31) 



-Plx + «S 



0, -wfj + -Pi x = 0, Pi = a fa. 
a 



We obtain solutions for the last equations in the form 

u±\x, y, z, t) = pi(x, y, z, t) + uf\y, z, t) 
Pi(x,y,z,t) =api(x,y,z,t) 



(32) 



where u{\y, z, t) is an arbitrary function. Without loss of generality we 
assume that Ui(y,z,t) = 0. 

Let us substitute ( |3~0~| in to (25)-(29). Equating expressions at e 2 in (25), 
(26), (29) and at e 3 ^ 2 in (27), (28) we obtain the following system of equations 
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in the case of x = 0(e 



1/21 



Pit 



Plx + u£? + fog? + 1U f] + = 0, 



(1) 

It 



-«S 2 J + - p ltf 

a 
a 

P 2 = ap 2 + a x p{ + (3p lxx - Hp 



(i) 

5 



a 



2.v 



(33) 



0. 



Using (32) we find from (33) 



Pit 



(2) 
It 



(3) 
"is 



1 + 

<*Ply, 
IPlz 



a 



/3 „ 
2a 



5 



(2) . 7 (3) 



(34) 



We can write (j34J) as the equation 
Pit + I 1 



tt'i 

a: 



PlPlx + 7^— Plxxx 

2a 



x 
2a~ 



Pla 



5 2 „ 7 2 - 
+ ~2~ Pira + y Piz* 



(35) 



Equation (35 ) is the three-dimensional counterpart of the Korteweg-de Vries- 



Burgers equation for description of nonlinear waves in the liquid with gas 
bubbles. 

In the case of 7 = (i.e. two-dimensional waves) and x = (i.e. the 



liquid is ideal) (35) reduced to the Kadomtsev-Petviashvili equation [19] . At 
7 = (35) is the two-dimensional Korteweg-de Vries-Burgers equation 20 



From relations (19) we see that the coefficient at the dispersion term in 
(35) takes the form (3 /2a = (piRqC^) / (6nP L) at a = and one can see that 
the value of dispersion depends on the radius of bubbles, the polytrophic 
exponent and the pressure in the unperturbed state of the bubble-liquid 
mixture. 

In the same way we can obtain that dissipation of the nonlinear waves 
depends on the liquid viscosity. In our consideration we do not take into 
account heat-transfer between the liquid and the gas in bubbles. However 
we can consider the influence of the heat-transfer on dissipation of nonlinear 
waves using the effective viscosity [l}[2]- This viscosity consist of the two 
parts: the liquid viscosity and the "thermal viscosity". 

The "thermal viscosity" depends on the thermal conductivity of the gas, 
the volume gas content and the radius of bubbles in the unperturbed state, 
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the frequency of bubbles oscillations and the polytrophic exponent [1}|2]. For 
example, from the formula for effective viscosity [TJ|2) we can see that in the 
case of a liquid with large viscosity and small bubbles the liquid viscosity 
causes dissipation of nonlinear waves. Otherwise dissipation of nonlinear 
waves is caused by heat-transfer and we can use the formula for the effective 
viscosity. 

Replacing in ^ the liquid viscosity on the effective viscosity we can take 
into account the dissipation of nonlinear waves due to the heat-transfer. 



4 The Painleve test for the thee-dimensional 
nonlinear evolution equation. 

Using transformations 



Pi 



f3 , 2a 



2{a + a,i 



P 




11/3, 1 p , 

-y , z = —\ — z 
7 V 2a 



(36) 



5 V 2a 



from (35) we have the following equation (the primes are omitted) 



(u t + uu x + u xxx - X u xx ) x + -(u yy + u zz ) = 0, 



(37) 



where X = f • 

Let us neglect the viscosity of a liquid. In this case x — and (37) takes 
the form 



(lit + uux + u xxx ) x + ^ (u yy + u z 



0. 



(38) 



To investigate the integrability of (38) let us apply the WTC-test with 



the Kruskal ansatz 21 22 . We are looking for solution of (38) in the form 



u 



(x,y,z,t) = ® P J2 A 3&> A, = A^y.z.t), $ = x + ^(y, z, t), (39) 

3=0 



where ip(y,z,t) is an arbitrary function. 

The leading terms in (38 ) Substituting (p9|) into these 



terms and equating coefficients at the lowest order of x + ijj(y,z,t) to zero 
we have 



p=-2, u = -12 



(40) 
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Substituting the expression 

u = -12(x + 1>(y, z,t))~ 2 + Aj(x + ij>{y, z,t)) j - 2 



(41) 



into leading terms of (|38j) and equating coefficient at Aj to zero we have the 
following Fuchs indices 



ji = -1, h = 4, h = 5, j 4 = 6. 



(42) 



We recall that (38 ) has the Painleve property if functions ip(y, z, t), A 4 (y, z, t), 



A 5 (y,z,t), A 6 (y 1 z,t) in (39) are arbitrary. 



Substituting (39) into 



38) and consistently equating terms at various 



powers of x + ip(y, z, t) to zero we find 
A 1 (y,z,t) = 0, A 2 (y,z,t) = -^t~«+^), A 3 = i {^ yy + ^ zz ) . (43) 

Functions A^(y, z, t), A 5 (y, z, t) can be taken as arbitrary. However the func- 
tion A^(y, z, t) can be arbitrary if the following condition is satisfied 



0. 



(44) 



From (44) we see that the function ip(y } z, t) cannot be taken as arbitrary. 
Thus (38) does not have the Painleve property. 

In the same way we can also show that (37) does not have the Painleve 
property. 



Exact solutions of the three-dimensional non- 
linear evolution equation 



We have obtained that (37) is not integrable in the general case and the 



Cauchy problem for this equation cannot be solved by the inverse scatter- 
ing transform. However there are exact solutions of this equation. Let us 
construct some of them. 

Using the traveling wave variables u(x, y, z, t) = U(9), 6 = x + y + z — C t 



we have the nonlinear ordinary differential equation from (37) 



C 2 + Ci 9 + (1 - C ) U + \u 2 + U ee - xU t 



(45) 



In the general case (45) cannot be solved but at C\ = (45) is equivalent 



to the traveling wave reduction of the Korteweg-de Vries-Burgers equation. 
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This equation was studied many times 20 23 26 . As this takes place differ- 
ent methods can be applied for finding exact solutions 27 31 . The general 
solutions of (45) at C\ — was obtained in [32]. It takes the form 



,.C--l + f-ex P {f}p(c 3 - 7 i-ex P {f},0,C 1 ),(4a) 

where p is the Weierstrass elliptic function, C3, C4 are constants of integra- 
tion. There is correlation on Co in the form 



C r = l±^C 2+ 3 ^ (47) 




(1 2) A V 

where Ck is defined by (47) as well. It is worth to note that this solution 



was obtained for the two-dimensional Korteweg-de Vries-Burgers equation 
in (20) . 

We see that (48) describe the weak shock wave. The amplitude and 



steepness of this wave defined by parameter \ which depends on the liquid 
viscosity. Waves of type (48) in the one-dimensional case were obtained 
as solutions of Korteweg-de Vries-Burgers equation in [6 33 . Weak shock 
waves of type (48) were observed experimentally as well [2|[5 ,33 . The plot 
of solution (48) at x = and t = is demonstrated in Figure 1. 

Let us consider equation (37) at x = 0. In this case we neglect the 
liquid viscosity. As we show above equation (37) is not integrable at % = 0. 



However we can construct multi-soliton solutions for this equation using 
special chose of coordinates. Following the Hirota bilinear method we look 



for the one-soliton solution of (37) in the form 34,35 



d 2 

u = A —— In F 

ox 1 

where F is the following 

F = l+exp{k x x + ky y + k z z — Q t + (f> }, Q 



(49) 



(50) 



It is worth to note that Q is defined by dispersion relation for (37). 



Substituting (49) and (50) into (37) we find that A = 12. Thus one-soliton 



solution of (37) has the following form 



u 



d 2 ( r i 

12 — In (1 + exp <yk x x + k y y + k z z - y 



t + 



(51) 



As we well know weakly nonlinear waves in bubbly liquids can be de- 



scribed by the Korteweg-de Vries equation 4 6,36 . We can consider solu- 



tion (51) as three-dimensional counterpart for the one-soliton solution of the 



Korteweg-de Vries equation. Further we shall show that the one-dimensional 
one-soliton solution is stable to transverse perturbations. 



We can find the two-soliton solution of (37) using formalism (49) as well. 



In this case the function F takes the following form 34 , 35 



F = 1 + exp{£;« x + y + z-Sl x t + 4> { q ) }+ 



+ exp{A;( 2 ) x + kf y + kf z-VL 2 t + <pf} + exp { + k^) x+ (52) 
+(kM + fcf) y + (k^ + fcW) z - + fl a ) t + 



<2, + 4 
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where 




Substituting Q and into (|37|) we find that A = 12 and $' 2) is 
determined by the following relation 



<A (1 ' 2) 
e <P 



° ) STx ) \Tx K x ) V^x S x V ) 

(*«*?> - W)'] [Q (k^) 2 (k^) 2 (*« + k^) 2 - (54) 
_ fi-(%( 2 ) - i-Wi-CO^ 2 - (k^kW - A-WfrWY 



Substituting (52) and (54) into (|49j) at A = 12 we have the two-soliton 
solution of (37). This solution is the thee-dimensional counterpart for the 
two-soliton solution of the Korteweg-de Vries equation. 
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The plot of the two-soliton solution of (37) at x = at various moments 



of time is presented in Figure 2. From Figure 2 we can see that solitary waves 



described by (37) interact elastically. 



The stability of the solitary wave solution for (|35|) at 7 
in a two-dimensional case) was shown in 14 



and x = (i.e. 



The full analysis of a stability 
of the solitary wave solution for (35) in linear approach is given in 37 at 
7 = 0. It is shown, that this solution is stable at 5 > 0. Let us show, that 



the results of 37 can be extended on (35) at 7 7^ 0. 



Without loss of generality we can transform ( 35 ) at x 
form 

(Au t - Au x + Auu x + u xxx ) x + 3u yy + 3u zz = 



to the following 



(55) 



We see that (|55|) has the stationary solitary wave solution in the form 

(56) 



3 cosh 



Following 37 we consider small perturbations of the solitary wave solution 

(57) 



u = u^ + e xt+ik y y+ikzZ v(x) 



Substituting (57) into (|55j) and linearizing the result we have 
v xxxx - 4[(1 - m (0) ) v} xx + AXv x + 3 {hi + hi) v = 0. 



(5f 



The spectrum of (58) and it's solutions were investigated in 37 assuming 



that A is a complex number. It was shown that at real k^ + kj > solutions of 



( 58 ) are bounded. Hence for real k y and k z solitary wave solutions governed 



by the equation (35) at x = are stable. 



6 Conclusion 

Let us formulate the results of this work shortly. The one-dimensional non- 
linear evolution equations for the description of waves in a liquid containing 
gas bubbles are well known but for the best of our knowledge there are no 
nonlinear evolution equations for description of nonlinear waves in a liquid 
containing gas bubbles in the three-dimensional case. In this work' we have 
considered the three-dimensional nonlinear waves in a bubble-liquid mixture 
taking into account the viscosity of liquid. We have obtained that the waves 
can be described by the generalization of the Korteweg-de Vries-Burgers 
equation and the Kadomtsev-Petviashvily equation. 

We have checked this nonlinear evolution equation on the Painleve prop- 



erty using the WTC test with the Kruskal anzats. We have found that (37) 
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is not integrable in the general case. However there are some exact solutions 
of (37) for description of nonlinear waves in a liquid with gas bubbles. We 
have presented traveling wave solutions for this equation. In the case of non- 
viscous liquid the nonlinear evolution equation is nonintegrable as well. But 
we have obtained one-soliton and multi-soliton solutions for this equation. 
We have shown that solitary waves governed by (37) in the absence of the 
dissipative losses is stable in relation to transverse perturbations. 
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